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Abstract 



C^ , New physical insight into the correspondence between path integral concepts 

a. 

and the Schrodinger formulation is gained by the analysis of the effective clas- 



sical potential, that is defined within the Feynman path integral formulation 
of statistical mechanics. This potential is related to the quasi-static response 
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1^ \ of the equilibrium system to an external force. These findings allow for a com- 

prehensive formulation of dynamical approximations based on this potential. 
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The path integral formulation of statistical mechanics is a powerful method to study 
quantum many-body systems. An essential property of this approach is the mapping of 
a quantum system onto a classical model of harmonic ring polymers, whose equilibrium 
properties can be studied with high accuracy. |1|] However, dynamical properties at finite 
temperatures can not be derived with the same sort of rigor, as the solution of the real time 
path integrals involves complex valued functionals without a suitable sampling function for 
stochastic integration. 

The effective classical potential (ECP), that was introduced by Feynman to study 
systems in thermodynamic equilibrium, [|| has been a central quantity to derive two quan- 
tum dynamical approximations, the quantum transition state theory (QTST), 0] that aims 
at calculating rate constants of activated processes, and the centroid molecular dynamics 
(CMD), that aims at calculating real time correlation functions for quantum particles. 
The definition of the ECP is based on the concept of path centroid, that is the center of 
gravity of each of the ring polymers that represent the quantum particle. However, in spite 
of the extensive use of the path centroid in condensed matter and chemical physics stud- 
ies, its relation to a measurable physical observable and hence the physical meaning of the 
derived dynamical approximations remain largely unexplained. 

In this work, we aim at a deeper understanding of the equivalence between the path 
integral and the Schrodinger formulation by showing that the centroid coordinate is related 
to the quasi-static response of the system to an external force. The Schrodinger formulation 
provides new physical insight into the dynamical approximations based on the ECP. At zero 
temperature, the ECP is the mean energy of minimum energy wave packets (MEWP's, to 
be defined below) and the CMD is an approximate dynamics based on these wave packets. 
Technical details of the analysis are left for a later work. [Q 

We begin by reviewing the definition of the ECP in the path integral formulation. The 
simple case of a quantum particle of mass m having bound states in a one- dimensional po- 
tential V{x) is considered. The Hamiltonian of the particle is Hq. The extension to the 
many-particle case is, for distinguishable particles, straightforward. |P] At a given temper- 



ature the equilibrium properties are derived from the partition function, Zq, and from the 
particle's probability density, p{x). Both are related by the expression Zq = J^^dxp{x). 
The path integral formulation of p{x) is 

r.=.m / S[X{U)] \ 

p{x) = D[x{u)] exp — , (1) 

Jx=x{0) \ "- / 

where S'[x('u)] is the functional of the Euclidean action of the path x{u), u is the imaginary 
time that varies between and (3%, and (3 is the inverse temperature {kBT)~^. The paths 
x{u) can be considered as alternative ways for the propagation of the particle, and the sum 
over paths can be analyzed in many different ways depending on the different classes into 
which the alternatives can be divided. The way conducting to the ECP uses the centroid 
or average point, Xc, of the path x{u) 

^c = ^^ du x{u) . (2) 

jjh Jo 

A class of paths is the subset of paths that have the same centroid. A constrained path 
integral over the class of paths with centroid at X is defined by introducing a delta function 
in the integrand of Eq. (|^), 

p^(x)=£D[x{uMX~x,)exp(-^^^] . (3) 

Px{x) may be considered as a probability density around the centroid position X. As 
illustration, the normalized function Z{X)~^px{x) obtained by a Monte Carlo path integral 
simulation of a particle of mass ?7i = 16 au in a double- well potential, V^p(x) = |(x^ ~ 1)^, 
is represented in Fig. ^ The probability density is shown for four values of the centroid 
coordinate X. We also show the probability density of the ground state of the potential 
Vdp{x) — fx, for several values of the parameter /, that represents an external force acting on 
the particle. The identity between these ground state probability densities and Z{X)^^px{x) 
will be explained below. The normalization constant Z{X) is an important quantity 

/■oo 

Z{X)= dxpxix), (4) 



which has the physical meaning of a probabihty density for the class of paths with centroid 
at X and it is called the centroid density. [^ The partition function can be recovered by a 
sum over all classes, i.e., 

/oo 
dXZ{X) . (5) 

-oo 

This integral, after substitution of Z{X) by the following definition, 

Z{X) = (^) exp [-f3F^f{X)] , (6) 

has the same form as the partition function of a classical particle moving in the potential, 
Fef{X), which is the ECP. We summarize some properties of the ECP: [^0 i) all the 
thermodynamic properties that depend on Zq can be derived from it; ii) it is temperature 
dependent; Hi) its calculation is analytical only for quadratic potentials, but variational 
approximations are available; iv) its high temperature limit is the actual potential V{x). 

The centroid density Z{X) (or the ECP) and the probability densities Z{X)^^px{x) 
are important quantities in the theory of path integrals, whose correspondence to the 
Schrodinger formulation has not been clearly stated. In the following, we show how the 
centroid density is related to a physical observable. We first consider a Hamiltonian de- 
pending on an external force / acting on the particle, Hf = Hq — fx. The path integral 
representation of the partition function Zf, corresponding to the Hamiltonian Hf, can be 
expressed as [§| 

/oo 
dXZ{X)e^f'' = Zo(e^^^) , (7) 

-oo 

where Z{X) is the centroid density for the particle with Hamiltonian Ho. The angle brackets 
show an average over the normalized centroid density. If Z{X) is known, the partition 
function Zf can be derived by Eq. (^ for any arbitrary value of the external force /. Thus, 
the analysis of classes of paths with fixed centroid allows to derive the thermodynamic 
properties of a whole family of systems whose Hamiltonian depends on a parameter /. The 
moments of the centroid density are defined as 



(X") = Z^^ / c/XZ(X)X" . (8) 

The physical meaning of Eq. (|^) is that the ratio of partition functions Zj/Zq is the 
function generating the moments of the centroid density. The moments are derived by 
differentiation of Zj/Zq with respect to the variable (3f. The result is simpler if we use the 
free energy, F/, defined as exp(— /5F/) = Zj. The first two moments are obtained as 



(^') - {Xf = -ksT {^\ . (10) 



9f' / /^o 



The mean-squared deviation of X is often called the classical delocalization of the particle. 
Higher order derivatives of Ff lead to higher moments of the centroid density. [Q If the 
moments are known, then the centroid density itself will be fully determined. We have then 
a formal relation between the centroid density and the change in the free energy as a function 
of a quasi-static force acting on the particle. This central result reveals the physical meaning 
of the centroid density within the Schrodinger formulation. 

This correspondence between the centroid density and the Schrodinger formulation is 
valid for arbitrary temperatures. The zero temperature limit is particularly interesting, as 
quantum effects are then most important. Before analyzing this limit, it is convenient to 
study a property of the Hamiltonian Hj that leads to the definition of MEWP's. We look for 
quantum states of the particle whose mean energy, E^^^ = (ipflHolipf), is minimum against 
small variations of {ipf). Moreover, the states must satisfy two constraints: i) their mean 
position is fixed at an arbitrary value x = {ipf\x\il!f); ii) they are normalized {'ipfl'ipf) = 1. 
By straightforward application of calculus of variations, Q one finds that \tpf) is the ground 
state of Hf 

{H,-fx)\ijf) = Ef\^f), (11) 

/ and Ef are Lagrange multipliers chosen so that the constraints i and ii are satisfied. Note 



that / is an implicit function of the arbitrary position x, i.e., / = /(x). The minimuni 
energy, ^J™"', as a function of x, is derived from Eq. (|1T]) as 

Er{x)=Ef + fx. (12) 

We call the states \ipf) the MEWP's of the potential V{x) whose mean energy is i?™"(x). 
We show next how these states are related to the ECP at zero temperature. 

At T = 0, the free energy is equal to the ground state energy, Ff = Ef, and from Eqs. 
(I), and ([10| ) we obtain 

W = -(^) ■,{X^)-{Xr = 0. (13) 

(X) is then the expectation value of the derivative of the Hamiltonian operator in Eq. (pUD 
with respect to / at / = 0, i.e., (X) is equal to the mean position of the ground state of the 
Hamiltonian Hq 

{X) = {iJo\x\tPo)=Xo. (14) 

The mean-squared deviation of the centroid density vanishes [Eq. (|T3|, thus the centroid 
density must be a delta function centered at xq. This result implies that only the class of 
paths with centroid at X = xq contributes to the path integral in Eq. (|lD, i.e., 

lim^(7V(a;) = hmZ(X)"Vx(a;) , for X = xq . (15) 

This identity provides a clear physical interpretation of path integral quantities at T = 0. 
It implies that the centroid density and the partition function of the particle are related 
by: Z{X) = Zq S{X — xq). Moreover, the asymptotic behavior of Zq and Z{xo) as T ^ 
is given by exponential functions of the ground state energy, exp(— /JE'o)) and the ECP, 
exp[— /5Fe/(xo)], respectively. Then, Fe/(xo) must be equal to the ground state energy. 
This result is a particular case of a relation valid for any arbitrary centroid position X. 
An essential step for this generalization is that the probability densities Z{X)''^px{x) are 
invariant under any change in the Hamiltonian Ho that is linear in the coordinate x. The 



final result of this generalization is: i) at T = the normalized probability density for 
a given class of paths, Z{X)~^px{x), is identical to the probability density of the MEWP, 
|(a;|^/)P, whose mean position is at x ^ X; ii) at T = the ECP, Fef{X), is equal to the 
mean energy, EJ^^^{x) of the MEWP whose mean position is at x ^ X. 

In Fig. |I| we showed the probability density Z{X)^^px{x) obtained from Monte Carlo 
path integral simulations of a particle in a double-well potential at temperature ksT = 10~^ 
au. This value is a small fraction of the lowest excitation energy, AEq = 29 x 10^^ au 
(tunnel splitting), and therefore a good approximation to the zero temperature limit. For 
comparison, the probability density of MEWP's, obtained numerically as the ground state 
of the potential V^p(x) — fx, are shown by broken lines. The identity of both probability 
densities is an important result that clarifies the physical meaning of fixed centroid path 
integrals. In Fig. |^ the energies Vdp{x) and £'™"(x) are shown. These curves correspond to 
the limits of infinite and zero temperature for the ECP, respectively. 

In QTST the second derivative of the ECP, Ff.f{X) with respect to X, has been shown 
to be an important quantity to determine the pre-exponential factor of rate constants. |Q If 
the effect of a small external force on a stationary quantum state is approximated by a rigid 
spatial displacement of the state, it can be shown [Q that in the zero temperature limit the 
second derivative of the ECP with respect to X gives an approximation to the first excitation 
energy of the Hamiltonian Hf, i.e. 
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'92E™"(x) 
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AEf . (16) 



m \ dx"^ 

In Fig. ^ we display the exact values of AEf as a function of x (remember that / is 
an implicit function of x). The l.h.s of Eq. ([16|) , drawn by a broken line, was obtained 
by numerical differentiation of the function EJJ^^^{x), shown in Fig. ^ At x = the 
approximation overestimates the energy of the tunnel splitting, AEq, by 25 %. 

An approximation to the time evolution of the particle can be formulated with MEWP's: 
i) only MEWP's with mean position x and momentum p are allowed as time dependent 
states, i.e., these are of the form (xlipf) exp{ipx/h). The value of the mean force is / = — /; 



a) the time dependence of x and p is given by the Ehrenfest relations 

dt m ' dt 

The total energy of the MEWP, £'™"(a;) +p^/(2m), is conserved along this time evolution. 
This approximation, formulated without any reference to path integral concepts, can be rec- 
ognized as the zero temperature limit of CMD. The formulation of CMD as an approximate 
dynamics based on MEWP's provides new physical insight into this approximation. In Fig. 
§ we show a phase space representation of the time evolution of x{t) and p{t) for an initial 
MEWP with x(0) = —0.7 au and p{0) = moving in the potential Vdp. The exact trajectory 
was obtained by numerical solution of the time dependent Schodinger equation. The CMD 
trajectory, derived from Eq. {^), is that of a classical particle moving in the ECP shown 
in Fig. ^. For comparison the phase space trajectory of a classical particle in the potential 
Vdp is also shown. The CMD result resembles the real time trajectory. 

Summarizing, the centroid density has been derived in the Schrodinger formulation with 
the help of a function generating its moments. The centroid density is related to the response 
of the system to a quasi-static force. For potentials with bound states the centroid density 
converges to a delta function as T ^ 0. Several results have been found in this limit: i) the 
ECP is the mean energy of the MEWP's of the potential; ii) the second derivative of the 
ECP with respect to the centroid position approximates the first excitation energy of the 
system; in) CMD is an approximate dynamics based on MEWP's. 

This work was supported by DGICYT (Spain) under contract PB96-0874. We thank E. 
Artacho and J.J. Saenz for helpful discussions. 
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FIGURES 
FIG. 1. The continuous lines show normalized probability densities, Z{X)^^px{x), obtained 
by fixed centroid path integral simulations at ksT = 10~^ au, in the potential Vrfp(x). The dotted 
lines (on top of the continuous ones) are the probability densities of the ground state of the potential 
Vdpix) — fx, that is shown by broken lines for several values of /. 

FIG. 2. The ECP in the limit T ^> (broken line) and the double well potential Vdp (contin- 
uous line). 

FIG. 3. First excitation energy ^Ef, for the Hamiltonian Hf (continuous line) as a function 
of the mean position x of the ground state of Hf, and its approximation based on the second 
derivative of the ECP (broken line) . 

FIG. 4. Phase space trajectory showing the time evolution of the mean position and momen- 
tum of an initial state corresponding to a MEWP of the Vdp potential with x{0) = —0.7 au and 
p{0) = 0. The exact result is compared to the CMD approximation and to the classical trajectory. 



10 





1.5 




-- 




0.5 







>^ 




^— > 


1 


•i-H 




C/3 




a 






0.5 


>. 









^B^ 




•i-H 




X) 




cd 




X) 




o 


0.5 


^ 




Oh 










0.5 - 







X = -0.71 au/ 

0.03 au / 
/ 



= - 0.41 au / 
0.01 au / 
I 



X = 0au 
f = Oau 



0.6 



- 0.3 








- 0.3 








- 0.3 








- 0.3 







3 
ui 

!=l 
f— I 

-4— ' 

!=l 

o 

Oh 



-2 



-1 







1 



2 



position (au) 










, , ■ 


1 


1 ' i ' 




0.4 


- 


1 \ 




/ 1 
/ 1 
/ 1 


- 






I \ 
1 \ 
I \ 


/-^ 


/ / 
/ / 
/ / 




0.2 


- 


\ V / 


r \ 


V y / 


- 




- 






\' / 


_ 


n 




1 Viy 1 


1 


1 Viy 1 





-2 



-1 







1 



X (au) 






0.4 - 



0.2 - 








-1.5 1 0.5 0.5 1 1.5 



X (au) 



IOh 



1.5 



1 



0.5 - 







-0.5 - 



-1 



-1.5 



1 


J. 1 1 1 1 1 1 1 1 1 

♦ ♦ 
* * 


_ ; 


• • 
• • 
• • 

■ \ classical 


; 


CMD 


■ 


;//^" ^^^"^"^vN 


1 


^^\ exact /y 


". 


• • 

• • 

• * ^ 


1 


*!*' 



-1.5 -1 -0.5 0.5 



1 



1.5 



x(t) 



